Quantum Theory for Generation of Nonclassical Photon Pairs by a Medium with 

Coherent Atòmic Memory 
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We present a fully quantum mechanical treatment of recent experiments on creation of collective 
quantum memory and generation of non-classically correlated photon pairs from an atòmic ensemble 
via the protocol of Duan et al. [Nature 414, 413(2001)]. The temporal evolution of photon num- 
bers, photon statistics and cross-correlation between the Stokes and anti-Stokes fields is found by 
solving the equation of motion for atòmic spin-wave excitations. We consider a low-finesse cavity 
model with collectively enhanced signal-to-noise ratio, which remains still considerably large in the 
free-space limit. Our results describe analytically the dependence of quantum correlations on spin 
decoherence time and time-delay between the write and read làsers and reproduce the observed data 
very well including the generated pulse shapes, strong violation of Cauchy-Schwarz inequality and 
conditional generation of anti-Stokes single-photon pulse. The theory we developed may serve as a 
bàsic approach for quantum description of storage and retrieval of quantum information, especially 
when the statistical properties of non-classical pulses are studied. 
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I. INTRODUCTION 



For transfer of quantum states between different nodes 
of quantum networks, the obvious choice is to employ 
photons as the fast and robust carriers of quantum in- 
formation [1]. However, future development of quan- 
tum communication intimately depends on successful at- 
tempts to reduce strongly the photon losses at large dis- 
tances that limit the range of application of this tech- 
nique. One possible way to cover large distances is us- 
ing quantum repeaters [2] , which combino the teleporta- 
tion of entangled states as a channel of quantum infor- 
mation transfer (cntanglement swapping) [3,4] with lo- 
cal atòmic memòries for quantum information. Among 
the various schemes for physical implementation of quan- 
tum repeaters and, hence, for realization of sealable long- 
distance communication [5-7] , one of the most promising 
approach is the protocol of Duan, Lukin, Cirac and Zollcr 
(DLCZ) [7]. This proposal is a probabilistic schcmc that 
relies upon entanglement between distant atòmic ensem- 
bles, which is created via successful detection of single 
photons cmitted by initially indistinguishable sources. 

In this scheme, by utilizing simple linear optical oper- 
ations, a quantum repeater protocol with built-in atòmic 
memory and entanglement purification is implemented, 
which is robust against the realistic imperfections, such 
as spontaneous loss and coupling inefficiency, that en- 
ables one to overcome photon attenuation and to commu- 
nicate quantum states over arbitrary long distances with 
only polynomial costs. In contrast to proposed earlier 
more complex protocols [2,4,8], the DLCZ scheme is at- 
tainable for current experimental technologies. This sig- 
nificant advance has been reached by exploiting the col- 
lective enhancement of atom-light interactions provided 



by optically thick atòmic ensembles due to many àtoms 
constructive interference effects [9-13] resulting in collec- 
tively enhanced signal-to-noise ratio. 

The initial step toward realization of DLCZ protocol 
is generation of an entanglement between remotely lo- 
cated atòmic ensembles. To that end the atòmic ensem- 
bles with all àtoms are prepared initially in the ground 
state | 1 ) (Fig la) are illuminated by a short-pulse weak 
Raman-pumping làser (referred to as a write làser), that 
induces a Stokes-photon emission in the transition | 3 ) 
— >| 2 ). Due to coherent coupling of different àtoms to 
the Stokes pulse propagating collinearly with the pump- 
ing làser, a collective atòmic mode is created in the form 
of long-lived spin-wave excitations. As a result, a strong 
correlation between these atòmic and Stokes-light modes 
is produced, while other optical modes are weakly corre- 
lated with the collective atòmic state and contribute to 
noise. This significant enhancement of signal-to-noise ra- 
tio allows one to generate a pure entanglement between 
two distinct atòmic ensembles by interfering and detect- 
ing single Stokes photons emitted from them. The chal- 
lcnging task for this scheme is to demonstrate that the 
collective atòmic state can be separately measured and, 
hence, the quantum correlation between the atòmic and 
signal-Stokes modes can be revealed. This can be done 
by first converting the atòmic spin excitation into state 
of the single-mode idler (anti-Stokes) light by applying, 
after a controllable time delay (Fig.lb), a second làser 
beam (read làser) , and then detecting the idler light again 
through a single photon detector. 

Very recently this program has been successfully ac- 
complished in a number of experiments on generation 
of nonclassical photon pairs both in samples of cold 
àtoms [14-16] and in room-temperature atòmic vapor 
celis [17,18]. An atòmic vapor celi has been used also 
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in another experiment [19] demonstrating convincingly 
a strong correlation between the Stokes and anti-Stokes 
pulses in the regime of largc photon number. Meanwhile, 
a theoretical description of these processes is absent with 
the exception of some particular results obtained for only 
the Stokes-photon cmission [7] or on the basis of phe- 
nomenological approach [14]. Since in the DLCZ scheme 
we deal only with spontaneous emission, a fully quantum 
treatment is needed. 

In this paper, we develop the quantum theory of cre- 
ation of atòmic memory and generation of nonclassical 
photon pairs via the DLCZ protocol. In our approach, 
to exploit the collcctive enhancemcnt of the atom-light 
coupling, we consider the interaction between àtoms and 
forward-scattered Stokes and anti-Stokes photons in low- 
finesse ring cavity [20] with an independent photon emis- 
sion by individual àtoms. The correlation between Stokes 
and anti-Stokes photons, separated by time interval con- 
sidcrably larger than the excited atòmic state lifetimc, 
is established via temporal dynamics of the collective 
atòmic mode. The equation of motion for the latter is de- 
rived by adiabatic elimination of cavity modes in the bad- 
cavity limit. The probability of excitation by the write 
and read làsers is assumed to be small and is treated in 
a perturbative way. Our results describe analytically the 
evolution of quantum correlations, including their depen- 
dence on coherent atòmic memory lifetime and on time 
delay between the write and read làser pulses. 

For future technological developments, new sources for 
single-photons are needed, which produce onc photon on 
demand, at a specific time and not at random. A re- 
markable feature of the DLCZ protocol is the ability to 
produce not only single-photon pulses, but also any spe- 
cific multiphoton state on demand via conditional mea- 
surement on quantum systems of correlated photon pairs. 
Compared to other systems of conditional generation of 
single photons, such as atòmic cascades [21,22] and para- 
metric down-conversion [23], this scheme has an advan- 
tage of simultaneous control over both the photon num- 
ber and the spatio-temporal shape of generated pulses 
that has been recently demonstrated in [15,16,18]. Our 
theory reproduces the prominent features of observed re- 
sults including pulse shapes and conditional generation 
of anti-Stokes single-photon pulse, as well as strong vio- 
lation of the Cauchy-Schwarz inequality. We believe that 
the calculations with this simple model may directly be 
applied to realistic situation of free-space atòmic celis. 
This expectation is supported also by an important ob- 
servation [24] that the three-dimensional theory confirms 
the large enhancement of signal-to-noise ratio predicted 
by simple cavity-QED models [7,11,25]. 

The outline of the paper is as follows. In section II 
we derive the master equation for the collective atòmic 
mode by eliminating the degrees of freedom associated 
with the cavity fields. The mean numbers of the Stokes 
and anti-Stokes photons inside and emitted from cavity 
are found in section III. Here we compare our results 
for photon fluxes with the experimentally observed pulse 



shapes. In section IV we calculate the photon statis- 
tics and cross-correlation between the fields and show a 
strong violation of Cauchy-Schwarz inequality. In section 
V we obtain the probability of conditional generation of 
single anti-Stokes photon and show its dependence on 
spin-decoherence time. The conclusions are summarized 
in section VI. 



II. D YN AMICAL EQUATION FOR COLLECTIVE 
ATÒMIC MODE 

We consider a large ensemble of N àtoms with a level 
structure shown in Fig. 1 and assume that initially all the 
àtoms are prepared by optical pumping in the ground 
state | 1 ). The sample in the state (Fig.lb) interacts 
with a weak write pulse with the duration TV- The write 
làser beam acts on the transition | 1 ) — >| 3 ) with a 
large detuning and Rabi frequeney Çlw and induecs 
spontaneous emission of a Stokes photon in transition | 3 
) — >| 2 ) whilc flipping an atòmic spin into the second 
ground state | 2 ). Detection of a forward propagating 
Stokes photon projeets the state of atòmic sample onto 
the nonclassical collective state í'l with the symmetric 
distribution of the fiipped spin excitations (see below). 
The read pulse acting on the transition | 2 ) — »| 4 } with 
the Rabi frequeney Q, R and duration Tr is applied after 
time delay and converts the stored spin excitations 
into the anti-Stokes light. 

Different schemes can be used for this retrieval pro- 
cess. It is possible to convert the stored spin excitations 
into an anti-Stokes photon in the electromagnetically in- 
duced transpareney (EIT [26]) configuration similar to 
that exploited in previous experiments [12,13] for restor- 
ing the classical light pulse. This mechanism has been 
used in the most of the experiments with DLCZ scheme 
[14-16,18]. Another way to transfer cfhcicntly a quantum 
state from àtoms to light is off-resonant Raman configu- 
ration demonstrated experimentally in [17]. Bcllow, we 
consider for simplicity the Raman configuration for emis- 
sion of anti-Stokes photon. Nevertheless, our results de- 
scribe very well the experimental data obtained with EIT 
configuration [15,18]. 

The fields of write and read làsers propagating along 
the axis of pencil-shape atòmic ensemble in z direction 
are given by 

Ew,n{t) = E w , R fH 2 R {t) exp(ik WtR z - ÍLü W . R t), 

(1) 

wherc fw,n{t) are the temporal profiles of the pulses. 

The positive-frequeney parts of forward propagating 
Stokes and anti-Stokes cavity fields at freqüències w\ and 
üü2 are expressed, respectively, in terms of annihilation 
(creation) operators ai(a^) and 02(0^") as follows: 
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, 2-küj., 



E\ + \z,t) = -i{—^-) 1/2 ai cy^{ihz - iüJit), i = 1,2 



(2) 



where V is a quantization volumc, which is assumed to 
be equal to the volume of the medium. 

Owing to the large detunings Ai = W31 — uw and 
A2 = W42 — üjr of the write and read làsers from the 
respective transitions, we can adiabatically climinate the 
upper states | 3) and | 4) , which are different in general. 
Then, under the conditions (kw — k\)L< 1 and (ku — 
k 2 )L < 1 , where L is the length of the atòmic sample, the 
interaction Hamiltonian for total system in the rotating 
framc has the form 



H = ih 



N 

E [ G <t)°: 

i=l 



Í£at-F(t)a$a 2 



+ h.c. (3) 



In Eq.(3) the summation is taken over all àtoms, a a a —\ 
a) i (P | is the atòmic spin-flip operator in the basis of 
the two ground states | 1) and | 2) for the i-th atom 
and 



F(t)=g AS -r 



9AS — J R (t), 



(4) 
(5) 



where the Rabi freqüències of classical fields are 
HsiEw/h and ÇIr = ha2Er/?i and 



9s 



1/2 



/i 32 and 9 as 



2TTLÜ2 

V 



1/2 



M41 



(6) 



are the atom - quantized fields coupling constants, pij is 
the dipole matrix element of the transition | i) — > | j). 
Note that the Stark shifts of the upper levels| 3) and| 4) 
induced by the write and read làsers are included in the 
freqüències of Stokes and anti-Stokes photons, so that 
u>\ = lüw —0J21 and lü2 = üjr +lü21- 

In terms of the collective spin 
operators 

/ 1 \ N 
the Hamiltonian (3) is written as 



H = ih^N [G{t)S + a+ - F(t)S+a 2 ] + h.c. 



(8) 



When all àtoms are prepared initially in level | 1), the 
states that are coupled by Hamiltonian (8) are totally 
symmetric. Particularly, if the atòmic ensemble is ini- 
tially in the state ^0 —\ li, I2, ....In), then upon emit- 
ting one Stokes photon it settles down into symmetric 

N 

statc * 1 = ^ £ I !i> ..2i, ....1 N ) (sce Fig.l). 



The system evolution is described by the master equa- 
tion for the whole density matrix p for the àtoms and 
cavity modes [27] 



d P 1 \tt 1 , ( d P\ 



(9) 



where the second term in right hand side (rhs) accounts 
for all relaxations in the system. With usc of the Liouville 
operator L[Ò]p = 2ÒpÓ+ - (Ò+Óp + pÓ+Ò) it is written 
in the form 



í d P\ 1 t\ 1 , T w(t) 

i=l i=l 



+ 



r«(t) 



N 



N 



E L ^2 ) ]p-7cE í(<) ^-/ , (°)]· 



i=l 



i=l 



(10) 



The first term in rhs of this equation, represents the cav- 
ity output at the freqüències Wj( i = 1, 2) with 2ki being 
the photon number damping rate for i-th mode, whilc 
the second and third terms describe spontaneous decay 
from the upper atòmic states |3) and| 4) into the ground 
levels |2) and |1) resulting in optical pumping (OP) 
to these states during the interaction with the write and 
read làsers, respectively. The rates of OP are 



IV(Í) 



-7^fw(t)l32, 



r«(t) 



A « (11) 



where 7^ is a partial decay rate of upper level i to the 
state j. 

We do not consider here two other channels (| 3 ) — >| 1) 
and| 4 ) — >| 2)) of spontaneous emission, because they do 
not change the atòmic spin distribution. The last term in 
rhs of Eq.(10) introduces into the model the relaxation of 
the atòmic ground-state coherence at the rate 7 C , which 
is supposed to be much smaller compared to the optical 
coherence damping rate 7. In this term jM = CT W _|_ 
represents the unit matrix in the basis of atòmic states. 

Upon introducing the operators of the atòmic popula- 

N 

tions of the ground states | 1) and | 2), S a =£ aiL ■, i — 

i=i 

1,2, we obtain from Eq.(10) the following equation for 
average number of the àtoms in the state | 2) : 



d(S 2 ) 



di { ±+2k)(nt ) -4 n) )+T 1 (t)N 



_[ 7c + r 2 (í)]<s 2 ) 



(12) 



where n. 



(in) 



1,2, are the mean pho- 



i=l 



ton numbers in, respectively, the Stokes and anti-Stokes 
modes inside the cavity. To simplify, the decay constants 
of photon numbers, k\ and fc 2 , are taken both equal to 
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k. Since the average number of output photons nf ut ^ is 
dctcrmined by the equations 



dn 



(out) 



dt 



= 2kn 



(*™) ■ _ 



we have 



d(S 2 ) 
dt 



dt [ni 



i = 1,2, 



4 íot) ) + ri(í)iv 



(13) 



-[ 7c + r 2 (í)] (s 2 ) 



(14) 



(tot) (in) 

n) ' — n) 



+n. 



(out) 



is the total number of 



where 
photons in i-th mode. 

The physical meaning of this equation is obvious show- 
ing that the atòmic population in the state | 2) is pro- 
portional to the photon number diffcrcncc in Stokes and 
anti-Stokes modes. Also, it increases due to OP from 
the state | 1), when the atòmic sample interacts with the 
write pulse, and this is described by the term TiN. Later 
on, it decreases at the rate T 2 proportional to the read 
pulse intensity, as well as due to escape of the àtoms from 
the làser beams area at the rate -f c . 

We now eliminate the cavity fields adiabatically treat- 
ing the two-photon interaction tcrms G(t) and F(t) in 
the Hamiltonian (8) in a perturbative way. By adiabatic 
elimination we obtain the equation of motion for the re- 
duced density matrix of the àtoms p a — Tr c p in the form 



= —L[S } Pa + —L[S] Pa + ( — ) rel , 



(15) 



where 



1 dt )rel ~ 2 



E 



+ 



r«(í) 



N 



N 



7< 



and 



2N 

a(t) = —G 2 (t) = af w {t), 
o/V 

^) = y F 2 (í)=/3/ fl (í) 



(16) 

(17) 
(18) 



are the Stokes gain and anti-Stokes absorption (from the 
state | 1)) coemcicnts, respectively. 

This is the central equation of our paper. In deriving 
Eq.(15), we have assumed only, apart from the weak in- 
teraction condition, that the write and read pulses are 
not overlapped in time. Using Eq.(15), we study in next 
scctions the quantum dynamics of the system including 
the evolution of collective atòmic mode, photon statistics 
and nonclassical correlations between the photons. 



III. MEAN PHOTON NUMBERS. PULSE SHAPES. 

In this section we apply Eq.(15) to find the mean pho- 
ton numbers in forward-scattered Stokes and anti-Stokes 
modes inside the cavity and the photon fluxes, given by 
Eq.(13), as a function of time, as well as the number of 
spin-wave excitations in the collective atòmic mode. 

For the overall atòmic population in the state | 2) we 
obtain from Eq.(15) 



d(S 2 ) 
dt 



a{t)[N sp (t) + 1] - f3{t)N sp (t) + Ti(t)N 



-[7c + r 2 (í)] (s 2 ) 



(19) 



where N sp ~ (S + S) is a number of spin-wave excitations 
or in other words, a number of the àtoms with flipped 
spin, which form the collective atòmic mode. It obeys 
the equation 



dN s , 
dt 



= a{t)[N sp (t) + 1] - (3{t)N sp (t) + Ti(í) - r tot 



(20) 



where T tot = lc + T 1 {t) +T 2 (t). 

In obtaining Eqs. (19,20), the commutation relation 



([S,S+]) cl 



(21) 



has been used, which follows from weak interaction condi- 
tion implying that upon interacting with the write làser, 
almost all the àtoms are maintained in the ground state 
|1>- 

The only diffcrcnce between the Eqs. (19,20) is due to 
the relaxation terms. Conscquently, the total number 
of the àtoms (S 2 ) in the state | 2) is essentially larger 
than the number of spin-wave excitations, because the 
former is basically determined by the optical pumping 
TiN from the ground state | 1), whereas the collective 
atòmic mode is generated at the rate a(t) « TiN as a 
result of coherent interaction with the write làser. At the 
same time, owing to the fact that the spin coherence and, 
hence N sp , decays at the OP rate of individual atom, 
the signal-to-noise ratio in Eq.(20) given by a/Vi and 
f3/T 2 is greatly enhanced due to the large factor of the 
atom number N [7] (see also bellow). 

From the comparison of the Eqs. (14) and (19) in the 
time intervals < t < T w and T 2 <t<T 2 +T Rl where 
T 2 = Tw + Td is the instance of read pulse switching on, 
we immediately find the equations for the mean photon 
numbers in cavity modes 



dn\ 



(in) 



dt 



= a(t)[N sp (t) + l]-2kn { ï 



(in) 



dn. 



(in) 



dt 



P(t)N sp (t) - 2kn 



(in) 



(22) 
(23) 



It is useful at this point to consider numcrical estima- 
tions. Reasonable parameters are: light wavelength A = 



800nm, fi 



w 



7, fl R ~ 10 7 , 7 ~ 10V 



7c 



10 4 s- 
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L = 1h-10cto, V ~ lcm 3 , and JV ~ 10 12 . In the free-space 
limit, k = c/L ~ 3 • 10 9 s _1 is the inverse of the propaga- 
tion time of the pulses through the atòmic sample. Then, 
from Eqs. (17,18) one has a ~ lQ 6 s~ 1 , (i ~ ÍO^" 1 , and 

ri ~ io- 5 7 , r 2 ~ io- 3 7 ~ 7c . 

Thus, the signal-to-noise ratio {a/T\ or /3/r 2 ) is about 
<~ 10 4 . At the same time a,(3 « k, so that from 
Eqs.(22,23) at kt » 1 we readily find 



n^ = ^a(t)[N sp (t) + l), 



4 n) = ^P{t)N sp {t). 
From the Eqs. (13) the photon fluxes are given by 



dn\ 



(out) 



dt 

dn { ° ut) 



dt 



a(t)[N sp (t) + 1], 



(24) 
(25) 

(26) 
(27) 



The solution of Eq.(20) with the initial value N sp (t — 
0) = has the form 



t t 
N sp (t) = J dt'a(t')cxp{ J dr[a(r) 



rtot(r)]} 

(28) 



It is worth noting that if one neglects the relaxations 
terms in Eq.(20), n^° ut \t) and N sp (t) obcy the same 
equation during the write pulse. This means that there 
is an unambiguous correspondence between the number 
of detected Stokes photons and spin excitations stored in 
collective atòmic mode. In particular, upon interacting 
with the write làser 



n<r t) {Tw)=N sp {T w ). 



(29) 



Simple expressions are found for rectangular làser pulses 



modes are the same ré{ wt \Tw) — n 2 out \T2 + Tr) ( we 
replace hereafter, for short, the time argument T 2 + Tr 
by T R ), provided that the intensity of read làser is suffi- 
ciently large, (3Tr >> 1, and the time delay between the 
làser pulses is shorter compared to the spin decoherence 
time: 7 C T 2 ~ j c T d « 1. 

In Fig.2 we show the Stokes and anti-Stokes pulse 
shapes calculated by means of the Eqs. (26,27) for the 
vàlues of parameters given above and for almost rect- 
angular làser pulses with the rise and fall times much 
shorter than the pulse duration. These results coincide 
with the theoretical calculations presented in [18] and 
reproduce very well the experimental data reported in 
[16,18], although in these experiments the anti-Stokes 
pulse has been retrieved in EIT configuration. Fig.2a 
demonstrates the transition from a spontancous to stimu- 
lated emission of Stokes photons with increasing of write 
pulse intensity. In Fig.2b, the flux of anti-Stokes pho- 
tons is depicted as a function of time for a fixed number 
of detected Stokes photons n^ ut \Tw) = 3. The total 
number of emitted anti-Stokes photons is determined by 
the areas of the corresponding peaks. It can be shown 
that in the case of strong read pulse, the stored spin 
excitations is completcly converted into the anti-Stokes 
photons, i.e. N sp {T w ) = n ( ° ut) (T w ) ~ n 2 ° ut) {T R ). This 
is evident also from Fig.3, wherc the evolution of spin- 
wave excitations is shown for different vàlues of read làser 
intensity with a fixed Raman scattering rate a(t) corre- 
sponding to n^ ut \T\y) = 3. Indeed, as it follows from 
Fig.3, after interaction with the read làser, N sp (Tr) = 
for ÇIr » Çlw, whereas a residual coherent excitation is 
preserved in atòmic ensemble, if Qr ~ Ctw- To demon- 
strate how the spin decoherence deteriorates the atòmic 
memory, we show in Fig.3b the same calculations for the 
case of ten times larger decoherence rate j c and with the 
same time delay between the write and read làser pulses. 
In this case the total number of retrieved spin excitations 
and, hence, of produced anti-Stokes photons is strongly 
reduecd. 



n[ rn \t) = ^-e at , n[ out \t) = e at - 1, < t < T, 
z/c 



and 



w 



(30) 



IV. PHOTON STATISTICS AND CORRELATIONS. 
VIOLATION OF CAUCHY-SCHWARZ 
INEQUALITY. 



( in )(i\ — _ft_„(° ut ) (rr„.\ z ,-/3(í-T 2 )-7cí 



(out)/,\ (out) /m \ 

n 2 (*) = n i ( T w) e 



fat 



: [l_ e -/3(t-T 2 )] ) T 2 <t<T 2 +T B 



(31) 



(32) 



where a and (3 arc defined in Eqs. (17) and (18). It is seen 
that always nf n) « nf ut) . From the Eqs. (30) and (32) 
it follows also that the numbers of output photons in both 



The nonclassical character of the Stokes and anti- 
Stokes fields generated in the DLCZ scheme has been 
experimentally studied in Refs. [14-18] by observing the 
violation of the Cauchy-Schwarz inequality. 

It is well known [28] that two electromagnètic fields, 
for which a positive true probability distribution exists, 
satisfy the following Cauchy-Schwarz inequality 



(12) 



<5 (11 V 22) , 



(33) 



where is the normalized second order auto- 

correlation functions for i-th field and g( 12 ^ is the cross- 
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correlation between the two fields. They are defined as 



,(«) 



(*) 



G(»)(t) _ ( at +(t)ot(t)oi(t)oi(t)) 
(n 4 (í)) 2 ~ ~ K(í)> 2 

G( 12 )(Í!,Í 2 ) 



i = 1,2 
(34) 



ff (12) (íl,*2) = 



(ni(íi)) (n 2 (í 2 )) 
(4fti)4(í 2 )a 2 ft 2 )ai(ti)) 
(ni(íi)) (n 2 (t 2 )> 



(35) 



with Uj(í) = (af (t)a,i(t)} , i = 1,2, being the mean 
photon numbers. 

The incquality (33) is violated for quantized fields. In 
our case, the correlation functions G^ (£) are easily cal- 
culated for the cavity modes. Then, taking into account 
that in the bad-cavity limit the output fields and cav- 
ity modes have obviously the same photon statistics, we 
apply the obtained results for the detected photons. 

Using the Hamiltonian (8), the Heisenberg-Langevin 
equations for the cavity modes are given by 



ai = VÏV^S++ J dt'F 1 (t')t 
o 



a 2 



y/N 



k 



S + 



z 

J dt'F 2 (t')e 



-fe(í-t') 



-k(t-t') 



(36) 



(37) 



where the noise operators Fi(t) associated with cavity 
fosses in the Stokes and anti-Stokes modes have the prop- 
erties [29] 



(fi(t)> = (Fimit')) = (F?mw) = o 



(38) 



{F i {t)F+{t'))=2k i 5 ii 5{t-t l ). 

The Eqs.(36) and (37) reproduce the solutions (24) and 
(25) for mcan photon numbers and n 2 n \ if the fol- 

lowing conditions: 

(S(t)F i (t')) = (F l + (t')S + (t)) = 0, i = 1,2 

(39) 

are satisficd for t > t' . By using the solution for S(t) ob- 
tained with the Hamiltonian (8) and applying again the 
properties of Langevin forces (38) it may be proved di- 
rectly that these correlations must vanish. Moreover, the 
similar calculations show that this is true also in general 
case of correlation functions with more than two field op- 
erators written in normal order. This allows us, keeping 
only the first terms in Eqs.(36) and (37), to express the 
correlation functions G(t) in terms of the atòmic collec- 
tive spin operators as 



a(t) 
2k 



*i(í), *i(í) = (S 2 (t)S+ 2 (t)) (40) 



G (22) W-(^) 2 $ 2 (Í), *,(*) = (S +2 (t)S 2 (t)) (41) 



The equations for $i(t) and <ï> 2 (í) are derived from the 
master cquation (15) 



rf$i(t < T w ) 
dt 

d® 2 (t > T 2 ) 
dt 



= [2a(í)- 7c ] $i(í), 
= -[2/3(í)+7c] * 2 (í), 



(42) 
(43) 



and have the following simple solutions 



$i(í) - $i(0) exp i j [2a{t') - 7c ]dt / , Í<T, 



(44) 



$ 2 (í) - *2(TV) exp i - / [2/3(0 + 7c ]dt' , Í>T 2 
i ^ I (45) 



Tw 



where the value $i(0) = (S 2 (0)S +2 (0)) = 2 is readily 
found recalling the commutation relation (21) and that 
N sp (0) = (5+(0)5(0)) = 0. To obtain ® 2 (T W ) we use 
the solution of the equation (28) at t = Tw 



Tw 



N sp (T w ) = (S+(T W )S(T W )) = cxp[ / a(z')dz'} - 1 

{ (46) 

which, together with <E>i(7V) from Eq.(44), yields 



$2(T W ) = (S +2 (T W )S 2 (T W )) = 2(S+(T W )S(T W )) 2 

(47) 

This relation indicates a chaotic nature of spin-wave 
bosonic excitations just like a similar relation (a +2 a 2 ) = 

2 (a + a) 2 takes place for a thermal light. 

By substituting the Eqs. (40,41, 44,45,47) and (24,25) 
for mcan photon numbers n^ n \Tw) and n 2 n \t) into 
Eqs. (34), we eventually have 



9 {11) (t <T W ) = 2 and ,g( 22 )(í > T 2 ) = 2exp[ 7c (í - T w )\ 

(48) 

showing that both the Stokes and anti-Stokes modes sat- 
isfy Gaussian statistics. It is worth noting that the corre- 
lations are independent of time apart from the last 
factor in g^ 22 ^ (i) indicating that an increase of time de- 
lay between the làser pulses or of spin-decoherence rate 
7c results in a superchaotic statistics of the anti-Stokes 
field. 

Similarly the cross-corrclation function G^ 12 \ti,t 2 ) 
may be represented in the form 



G^\t ut2 ) = ^^ 12{ tut 2 y, 

*i 2 (ti,*2) - (S(íi)S+(í 2 )S(í 2 )S+(íi)) 



(49) 
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where t\ < TV and t 2 > T 2 . The equation for $12(^1, t 2 ) 
is dcduced from Eq.(20) by applying the quantum regres- 
sion thcorcm [27,30] with t 2 = íi + t, where r > t^, 
yiclding 



= -[/?(íi+r)+7 c ] 0i 2 (ti,ti +r); 



(50) 



the solution of which is 



*1+T 



*i2(íi,íi+r) = 0i 2 (íi,ti)exp[- / 0{t')dt' - -y c r] 

íi (51) 

Using Eq.(47), $i 2 (íi,íi) is easily transformed to 



*i 2 (ti,íi) - (S^S+ih)) [2(S+(í 1 )S , (ti)) + l] 

(52) 

From Eq.(33), the cross-correlation between the two 
modes is then obtained as 



V. PRODUCING SINGLE PHOTONS VIA 
CONDITIONAL MEASUREMENT. 

The strong non-classical correlation between the 
Stokes and anti-Stokes pulses in DLCZ scheme can be 
employed for producing single photons as a result of con- 
ditional measurement on correlated photon pairs (see also 
[22,23,31]), when a detection of one Stokes photon tightly 
projects the anti-Stokes field into an one-photon state. In 
order to quantify the confidence level of this procedure 
we analyze the third-order correlation function [28] 



i (3) (h,t 2 ,t 2 ) 



Pc(h,h I li ) 

1^(12 I ll)| 



2 ' 



(55) 



which is less unity for quantized fields (it approachcs 
zero in the case of creation of single photons), while 
g^\ti,t 2 ,t 2 ) > 1 for classical fields. In Eq.(55), 
-Pc(l2,l2 | li) and P c (l2 | li) arc conditional probabili- 
ties for detection of two and one photon from anti-Stokes 
field, respectively, conditioned upon the detection of one 
photon in the Stokes channel. Using the known formulac 
[28,30], the Eq.(55) may be transformed to 



,(12) 



(Í1.Í2) 



(S + (h)S(h)) 



(53) 



and, thus, for the ratio between |<?( 12 )| 2 and g^ 11 ^ ■ í^ 22 -* at 
íi = Tw and t 2 = T 2 + Tr we fmaily have 



7 (12) 



R = 



{ti,t 2 )\ 



l + 2n[° ut) (T w ) 



2n[° ut, {T w ) 



exp(- 



IcTd) 

(54) 



where (S + (Tw)S(Tw)) is replaced by the number 
rii * (TV) of Stokes photons emitted from cavity (see 
Eq.(29)). Remind that in deriving Eqs.(48) and (54) we 
have neglected all atòmic relaxations during the interac- 
tion with làser pulses under assumptions T tot ««,,(?, but 
we have kept the term j c Td, which is by no means small 
and leads to a loss of atòmic memory that is manifcsted 
by exponential decay of the ratio R. 

It follows from Eq.(54) that if the writc pulse is 
sufhciently weak, so that the excitation probability p 
in the Stokes channel is less than unity and, hence, 
n < f ut \Tw) — aTw = P <s£ 1, the Cauchy-Schwarz in- 
equality is strongly violated by the law R ~ 1/p 2 3> 1. 
In the experiments [14-17] the violation of this incquality 
has been tested to confirm the non-classical correlations 
between the Stokes and anti-Stokes pulses and a signifi- 
cant violation has been observed for small p, although, as 
has been discussed in [14,15], several impcrfections limit, 
in practice, the degree of violation of Cauchy-Schwarz 
inequality and lead to deviation from the ideal case ex- 
pressed by Eq.(54). 



P(ll)p(ll,l 2 ,ll) 



i^(ii,i 2 )r 

where the total probabilities are given by 

P{l 1 ) = (a+{t)a 1 {t))=n { r\t l ) 



(56) 



P(li,l 2 ) = G (12) (íi,Í2) 



a{ti)P(t 2 ) 
2k 2k 



x (S(h)S + (h)) exp[- / p{t')dt' - 7c (í 2 - íij] 

l (57) 

P(li,l 2 ,l 2 ) = G( 122 )(í 1 ,Í2,í 2 ) 

= (af (íi)o 2 "(í 2 )aj' (Í2)a 2 (Í2)a 2 (í 2 )ai(íi)) 

Here we have kept only the second term in G < - 12 ^(íi, t 2 ) 
(see Eqs.(49) and (52)) taking into account the smallncss 

of n^Hh) - (S+(íi)S(ti)). 

The last correlation function in Eqs.(57) is calculated 
by the same method employed in previous sections 

^ 22 )( Í1 , Í 2, Í 2) = ^(^) 2 ^ 1 , Í 2,Í2), 

*(ti, í 2 , t 2 ) - (S(t 1 )S +2 (t 2 )S 2 (t 2 )S + (t 1 )) , (58) 

where t\ < T^y, t 2 > T 2 . We obtain the equation for 
Í2, t 2 ) by applying again the quantum regression 
theorem to the Eq.(43) and find the following solution 
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Í2 



*(ti, Í2, Í2) - ti, íi) exp | - j [2/3(t') + 7c ]dt' 



where, for $(<i,ii,ii), by solving the corresponding 
equation obtained from Eq.(15), we have 



$(íi,íi,íi)=4[r4 ouí) (íi) + l] 



3 (out) I , \ H 



(59) 



Substituting Eqs.(57) and (58) into Eq.(56), we finally 
gct 



gW(T w ,T R ,T R ) = 4n[° ut \T w ) 

3 (out) 



exp(7 c r d ) (60) 



In the regime of very weak Raman excitation, 
n (° ut ) (x w ) ~ p <c 1, and providcd that j c Td < 1 one 
has g^ ~ 4p <C 1 that corresponds to ideal case of sin- 
gle photon crcation. 

Expcrimentally, the conditional generation of sin- 
gle photons in the anti-Stokes pulse has been recently 
demonstrated with g^ ~ 0.3 in an ensemble of cold Cs 
àtoms [15]. The possibility of conditional preparation of 
Fock states with a given number of photons in optically 
dence medium of Rb àtoms has been reported in Ref. [18] . 



VI. SUMMARY 

In summary, we have developed a theory in the weak 
pumping limit to describe analytically the non-classical 
correlations between the photons emitted separately from 
a medium with macroscopic atòmic memory for single- 
photon fields. First, the spin excitations in àtoms are 
produced in a Raman process induced by the write pulse 
and arc cntanglcd with forward-scattered Stokes photons, 
and, finally, they are converted by the read pulse into 
the anti-Stokes photons with an cfficicncy that depends 
on the read làser intensity. An important result of our 
study is that the correlations between the photons are 
purely determined by only the atòmic spin correlations. 
This allows one to describe the evolution of the system 
by the dynamical equation for the àtoms, instead of solv- 
ing the underlying equations for quantized pulses. We 
have derived the master equation for atòmic dynamics 
and have found its analytical solutions for consecutive 
write and read pulses with time separation longer than 
the duration of the pulses. The time dependence of pho- 
ton correlations shows that all quantum effects (violation 
of Cauchy-Schwarz inequality, conditional generation of 
single photons, etc) are robustly manifested within the 
atòmic memory lifetime. Otherwise, they exponentially 
disappear at the atòmic spin decohcrcncc ratc. To min- 
imize the dissipation we considered the Raman scheme 



for mapping the atòmic spin state onto the anti-Stokes 
light. Our results obtained in bad-cavity limit are consis- 
tent with those observed in single pass scheme with dense 
atòmic media, where for retrieving of the anti-Stokes light 
the EIT configuration has been used [15,16,18]. From 
this point of view, the theory we developed may serve 
as a bàsic approach for quantum description of storage 
and retrieval of quantum information, especially when 
the statistical properties of non-classical pulses are stud- 
icd. 

This work has been supported by the ISTC Grant No. 
A-1095. 
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FIGURE CAPTIONS 

Fig.l. a) Atòmic level structure for emission of Stokes 
and anti-Stokes photons in off-resonant Raman config- 
uration with one photon detunings Ai and A2, respec- 
tively. b) Interaction diagram in the case of one Stokes- 
photon emmission, where in initial state í'o a ll àtoms 
are prepared in the ground state | 1 > . In intermediate 
stage, the atòmic ensemble is in a collective state ^1 with 
symmetric distribution of one atom excitation. c) Time 
sequence for the write and read pulses having durations 
Tw and T#. is time-delay between the pulses. 

Fig.2. a) Stokes photon flux dji^ ut ^ /dt as a function 
of time for different vàlues of write pulse intensity. The 
curves are labeled with the corresponding vàlues of a. 
The total number of Stokes photons emitted from the 
cavity are determined by the areas of the corresponding 
peaks. The duration of write pulse T\y = 1-6/is. b) Anti- 
Stokes photon flux dn^ 11 ^ / dt for fixed number of output 

Stokes photons nf ut \Tw) — 3 and different vàlues of 
read làser intensity. The read pulse duration Tr = ífxs. 
The time-delay Td — 1.4/zs and j c — 0.03 (^s) -1 . 

Fig.3. a) Average number of spin-wave excitations 
N sp = (S + (t)S(t)) as a function of time for fixed number 

of detected Stokes photons nf ut \T\y) — 3 and the vàl- 
ues of read làser intensity, as in Fig.2b. j c = 0.03 (/xs) _1 . 
b) The same, but j c — 0.3 (^s) _1 . 
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